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^ , We present a numerical study of the structural correlations associated to gas adsorp- 

■ tion/desorption in silica aerogels in order to provide a theoretical interpretation of scattering exper- 

' iments. Following our earlier work, we use a coarse-grained lattice-gas description and determine 

the nonequilibrium behavior of the adsorbed gas within a local mean-field analysis. We focus on the 
^ ■ differences between the adsorption and desorption mechanisms and their signature in the fluid-fluid 

and gel-fiuid structure factors as a function of temperature. At low temperature, but still in the 
regime where the isotherms are continuous, we find that the adsorbed fluid density, during both 
00 . filling and draining, is correlated over distances that may be much larger than the gel correlation 

length. In particular, extended fractal correlations may occur during desorption, indicating the 
existence of a ramified cluster of vapor filled cavities. This also induces an important increase of 
the scattering intensity at small wave vectors. The similarity and differences with the scattering of 
fluids in other porous solids such as Vycor are discussed. 
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In a series of recent papers 0,0, 13 , we have presented a comprehensive theoretical study of gas adsorption/desorption 
in silica aerogels, revealing the microscopic mechanisms that underly the changes in the morphology of the hysteresis 
loops with temperature and gel porosity. In particular, we have shown that the traditional capillary condensation 
scenario based on the independent-pore modelQ does not apply to aerogels, as a consequence of the "open" nature and 
interconnectedness of their microstructure. We have found, on the other hand, that nonequilibrium phase transitions 
(that differ on adsorption and desorption) are at the origin of the very steep isotherms observed with ^He in high 
porosity gels at low temperaturej^ Ig]. In this work, we complete our study by investigating the correlations within 
■ the adsorbed fluid and computing the fluid-fluid and solid-fluid structure factors that can be measured (at least 
' indirectly) in scattering experiments. Scattering methods (using x-rays, neutrons and visible light) are now frequently 
1 combined with thermodynamic measurements for extracting information on the structure and the dynamics of the 
' adsorbed molecules and understanding the influence of solid microstructure on fluid properties,2j- In the case of 
''He in aerogel, both small-angle x-ray scattering (SAXS)0 and light scattering measurements^ have been recently 
performed along the sorption isotherms. However, the interpretation of the scattered intensity is always complicated 
by the fact that it contains several contributions that cannot be resolved without assuming some mechanism for the 
sorption process. For instance, in the case of a low porosity solid like Vycor, the evolution of the scattered intensity 
along the capillary rise is usually interpreted in the framework of an independent-pore model, with the gas condensing 
I : in pores of increasing size that are (almost) randomly distributed throughout the material [lOl fill . This explains 
that long-range correlations are not observed during adsorption. On the other hand, we have shown that large-scale 
Ch , collective events occur in aerogels, and this may have a significant influence on the scattering properties. Indeed, we 
O ■ shall see in the following that scattering curves at low temperature may not reflect the underlying microstructure of the 
gel. More generally, our main objective is to understand how the different mechanisms for adsorption and desorption 
reflect in the scattering properties as the temperature is c hang ed. (Note that there has been a recent theoretical study 
of this problem that is closely related to the present onep^: there are, however, significant differences that will be 
commented in due place.) Particular attention will be paid to the 'percolation invasion' regime that is predicted to 
occur during the draining of gels of medium porosity (e.g. 87% ) an d that manifests itself by the presence of fractal 
correlations. Such correlations have been observed in Vvcor [lol Ull fl^ and xeroeelfT^. but no experiment has been 
carried out so far to detect a similar phenomenon in aerogels. We therefore hope that the present work will be an 
incentive for such a study. On the other hand, the influence of gel porosity on scattering properties will only be evoked 
very briefly. In particular, for reasons that will be explained below, the correlations along the steep (and possibly 
discontinuous due to nonequilibrium phase transitions) isotherms observed in high porosity gels at low temperature 
are not investigated. 

The paper is organized as follows. In section 2, the model and the theory are briefly reviewed and the 
cart.texomputation of the correlation functions and the corresponding structure factors is detailed. The numeri- 
cal results are presented in section 3. The relevance of our results to existing and future scattering experiments is 



2 



discussed in section 4. 



II. MODEL AND CORRELATION FUNCTIONS 
A. Lattice-gas model 

As discussed in previous papersQ, 0, 0, our approach is based on a coarse-grained lattice-gas description which 
incorporates the essential physical ingredients of gcl-fluid systems. The model Hamiltonian is given by 

<ij> <ij> i 

where = 0, 1 is the fluid occupation variable (i — 1...N) and rji — 1,0 is the quenched random variable that 
describes the solid microstructure (when r]i = 0, site i is occupied by the gel; 4> — (1/-/V) Vi is thus the gel porosity). 
Specifically, we adress the case of base-catalyzed silica aerogels (typically used in helium experiments) whose structure 
is well accounted for by a diffusion limited cluster-cluster aggregation algorithm (DLCA)^^. In the Hamiltonian, 
Wff and Wgf denote respectively the fluid- fluid and gel- fluid attractive interactions, /i is the fluid chemical potential 
(fixed by an external reservoir), and the double summations run over all distinct pairs of nearest-neighbor (n.n.) sites. 

Fluid configurations along the sorption isotherms are computed using local mean-field theory (i.e. mean-field 
density functional theory), neglecting thermal fluctuations and activated processes (the interested reader is referred 
to Refs.P,|3,l3 for ^ detailed presentation of the theory). As fi varies, the system visits a sequence of metastable 
states which are local minima of the following grand-potential functional 

^{{Pi}) = ksTY^lp, In Pi + im - Pi) Mm - Pi)] - Wff Y P^Pi ~ ^af Y t^'^^ ~ '^j) + ~ ^ P^P' 

i <ij> <ij> i 

where PiUrji}) =< TiVji > is the thermally averaged fluid density at site i. Earlier work has shown that this approach 
reproduces qualitatively the main features of adsorption phenomena in disordered porous solids*16l Il7j. 

All calculations presented below were performed on a body-centered cubic lattice of linear size L = 100 {N — 2L^) 
with periodic boundary conditions in all directions (the lattice spacing a is taken as the unit length). L is large 
enough to faithfully describe gels with porosity (j) < 95%, but an average over a significant number of gel realizations 
is required to obtain a good description of the correlation functions. In the following, we use 500 realizations. Wff 
is taken as the energy unit and temperatures are expressed in the reduced unit T* — T/Tc, where Tc is the critical 
temperature of the pure fluid (Tc = 5.195K for helium and kTc/wff = 2 in the theory). The interaction ratio 
y = Wgf /wff is equal to 2 so as to reproduce approximately the height of the hysteresis loop measured with *He in a 
87% porosity aerogel at T = 2.42is:[3||]. 

B. Correlation functions and structure factors 

In a scattering experiment performed in conjunction with gas adsorption/desorption, one typically measures a 
(spherically averaged) scattered intensity I{q) which is proportional to a combination of the three partial structure 
factors Sgg{q),Sff{q) and Sgf{q), where g and / denote the gel and the fluid, respectively. Whereas a macroscopic 
sample is usually considered as isotropic and statistically homogeneous, our calculations are performed in finite samples 
and on a lattice, and some work is needed to obtain structure factors that can be possibly compared to experimental 
data. 

Let us first consider the gel structure factor. The calculation of Sgg{q) proceeds in several steps. As in the case 
of an off-lattice DLCA simulation [l^. we first compute the two-point correlation function gggiv) — hgg{r) -I- 1 by 
performing a double average over the lattice sites and the gel realizations, 

P^gggi^) = ^ X] (1 ^ '^'^^l ^ %)^r,r., - PgSr,0 (3) 

taking care of the periodic boundary conditions. Here, — — r^, pg = 1 — (j) is the lattice fraction occupied by 
the gel, and the second term in the right-hand side takes into account the fact that there are only one particle per 
site, which yields the (point) hard-core condition, ggg(r = 0) = 0. In this expression and in the following, the overbar 
denotes an average over different gel realizations produced by the DLCA algorithm. The computation of ggg (r) is 
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efSciently performed by introducing the Fourier transform Pg{q) = ~ ^j ) exp(— 2t7rq.ri /N) where q is a vector of 
the reciprocal lattice, and by using the fast Fourier transform (FFT) method|la.ri9|. This reduces the computational 
work to O(A^lnA^) instead of 0{N'^) when the direct real-space route is used. (The same method is applied to the 
other correlation functions.) 

In a second step, we "sphericalize" the correlation function by collecting the values having same modulus of the 
argument r, 

= (4) 

Finally, instead of storing the values of ggg{r) for all possible distances rij on the lattice between d = a\/3/2, the 
nearest-neighbor distance, and L/2, we bin the data with a spacing Ar = 0.05 and interpolate linearly between 
two successive points (the restriction to r < L/2 avoids boundary artefacts). Moreover, we impose the "extended" 
hard-core condition ggg(r) = for r < d, in line with our interpretation of a gel site as representing an impenetrable 
silica particle ^ll. (In Ref. 13^, in contrast, the model is thought of as a discretization of space into cells of the size 
of a fluid molecule and the gel particle "radius" is varied from 2 to 10 lattice spacings.) Of course, the interpolation 
procedure does not completely erase the dependence on the underlying lattice structure, especially at short distances. 
Following Ref. fl8| , the structure factor is then computed using 

sin(gr) 

Sgg{q) = 1 + 47rpg / r^ihggir) - hgg) dr (5) 

JO 1^ 

where hgg is a very small parameter adjusted such that Sgg{q) ^ as g ^ 0. Indeed, since the DLCA aerogels are 
built in the "canonical" ensemble with a fixed number of gel particles Ng = PgN, the following sum-rule holds: 

PsEffffsW^^s-l (6) 

r 

which readily yields Sgg{0) = 1 + PgJ2r^Jmi^) = in Fourier space. This trick allows one to obtain a reasonable 
continuation of Sgg{q) below q-min = 27r/L|18j. 

Similarly, the fluid- fluid two-point correlation function gff{r) ~ 1 + hff{r) is defined as 



where p/ = {'^i'Hi) ~ {^/^)J2i7'i the average fluid density and the sum is performed again over all 

lattice sites to improve the statistics (for notational simplicity, we have dropped the overbar on pf). Because of the 
double average over thermal fluctuations and over disorder, there are two distinct contributions to hff{v), which are 
usually called "connected" and "blocking" or "disconnected" [20lf2l|. and which, in the present case, are given by the 
expressions, 



N 

2, J 
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f^ffA^) = ^ Yj^PiPj - P}] ^r,r., • (9) 



In the pure fluid, /i//^c(r) is just the standard connected pair correlation function whereas hjf^dij) has no equivalent. 
It turns out, however, that only hff^di^^) can be computed along the sorption isotherms. Indeed, the quantity 
{TiViTjVj) cannot be obtained in the framework of mean-field theory, and the only available route to hff^c{i^) is via 
the "fluctuation" relationplj 

where /i^ is a site-dependent chemical potential|22j. However, this relation only holds at equilibrium (like the Gibbs 
adsorption equation pj = —{l/N)dVl/dp discussed in Ref.P|) and therefore it cannot be applied along the hysteresis 
loop where the system jumps from one metastable state to another. (In a finite sample, the grand potential changes 
discontinuously along the adsorption and desorption desorption branchesP|.) We are thus forced to approximate 
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hff(r) by its disconnected part, ^//,(i(r)|23]- However, this may not be a bad approximation at low temperature 
because the local fluid densities Pj are then very close to zero or one0, which likely implies that /i//^c(r) is a much 
smaller quantity than /i//,d(r)|2J|. 

We then apply to hff(r) the same procedure as for ggg(r), taking radial averages and then performing a binning 
of the data and a linear interpolation. There are no "extended" hard-core in this case. Indeed, since the scale of the 
coarse-graining is fixed by the size of a gel particle (typically, a few nanometers), a lattice cell may contain several 
hundreds of fluid molecules which may be thus considered as point particles. hff{r) is then also interpolated between 
r = and r = d. 

In a grand-canonical calculation, the number of fluid particles fluctuates from sample to sample, which implies the 
following sum-rule for the disconnected pair correlation function /ty_f^('r')|2ll| (and thus for hff{r) in our approxima- 
tion) 



Y^hMr)^T.h.ir)^N'-I^^^!i^ (11) 

where p/({?7i}) = (l/-^)X]iPi the average fluid density for a given gel realization. This sum-rule can be used to 
also extrapolate Sff{q) below q < 2n/L, using 

Sffiq) = l + 47rpf f'\\hff{r)~hjf)'^^dr (12) 
Jo 9' 



where hff is adjusted so that Sff{0) = 1 -|- Af[pj({?/i}) — Pf]/Pf- 

Finally, the gel- fluid two-point correlation function ggj (r) = 1 + hgf (r) is computed from 



PgPf9gf[^) = ^ 



W = 4E(l-^*)(^J-^^)'^r,r., (13) 



and then sphericalized, binned, and linearly interpolated (taking ggf(r) = for < r < d/2 since no fluid molecule 
can be found at a distance less than d/2 from the centre of a gel particle). The cross structure factor Sgf{q) is then 
obtained from 

sin for") 

Sgfiq) = An^p^pJ / r^{hgf{r) - hgf)-^dr (14) 
Jo 9' 

where hgf is adjusted so as to satisfy the sum-rule Sgf{q — > 0) = which again results from the absence of fluctuations 
in the number of gel particles. 



III. NUMERICAL RESULTS 



A. Aerogel structure 

We first concentrate on the case of the empty aerogel. We have already presented in Ref. 0] the pair correlation 
function ggg{r) for several porosities between 87% and 99%. These curves exhibit a shallow minimum that strongly 
depends on and whose position gives an estimate of the gel correlation lenght £,g, as suggested in Ref. 18]. A DLCA 
gel can indeed be sketched as a disordered packing of ramified blobs with average size £,g. For instance, £,g varies 
approximately from 4 to 10 lattice spacings as the porosity increases from 87% to 95% (this is much smaller than the 
box size L, which ensures that there are no finite-size artifacts in the calculation of the gel structure p5l|'). Only the 
highest-porosity samples exhibit a significant power-law regime ggg{r) ^ r~^^~'^f'> that reveals the fractal character 
of the intrablob correlations. Although we shall essentially focus in the following on the case of the 87% porosity gel, 
for the sake of completeness we show in Fig. 1 the evolution of the simulated gel structure factor with porosity. 

The curves closely resemble those obtained with the continuum model 18]. In particular, they exhibit the same 
damped oscillations at large q that result from the "extended" hard-core condition ggg (r) — Oioi r < d (the oscillations, 
however, are more significant in the continuum model). The range of the linear fractal regime increases with porosity 
(it is almost non-existent in the 87% gel) and corresponds asymptotically to a fractal dimension df « 1.9 (the value 
df ~ 1.87 was obtained in Ref.[l| from the ggg{r) plot for the 99% aerogel). A characteristic feature of the curves 
is the existence of a maximum at smaller wavevectors whose location q^ decreases with porosity and correlates well 
with 1/fg (q™ 2.6/^g). This maximum is thus the Fourier-space signature of the shallow minimum observed in 
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FIG. 1: Gel structure factors Sgg{q) obtained with the DLCA algorithm for different porosities. From left to right: <j!> = 
0.99, 0.98, 0.97, 0.95, 0.92, 0.90, 0.87. The dashed line has a slope -1.9. The arrow indicates the wavevector q = 2-k/L ^ 0.063. 
(Color on line). 



ggg{r). (Note that varying the size L has only a weak influence on the small q portion of the curves for the 87% and 
95% gels, which validates the continuation procedure used in Eq. 5.) 

To compute the scattering intensity /(g) and compare to the results of small-angle x-rays or neutron experiments, 
it is necessary to introduce a form factor F{q) for the gel particles. One can use for instance the form factor of 
spheres with radius R = d/2, F{q) = 3[sm{qR) — qRcoa{qR)]/{qR)^ . The curve I{q) = S^{q)F{q)'^ then differs from 
Sgg{q) in the large-g regime {q/2n > d~^) where it follows the Porod law /(g) ~ g~'*|2a|- other hand, the 

intermediate "fractal" regime (^^"^ < g/27r < d^^) where Sgg{q) ^ q^"^^ , and the location g^ of the maximum are 
essentially unchanged. By cornparing the value of g™ in Fig. 1 with the actual value in the experimental curves 
(displayed for instance in R,ef.|l8|). we can thus fix approximately the scale of our coarse-graining. For (j) = 95%, 
qm ~ O.OlA^^ which yields a « 3 nm, a reasonable value for base-catalysed silica aerogels which is in agreement with 
the estimation of Ref.P| obtained from the gel correlation length ^g. It is worth noting that the DLCA simulated 
structure factors present a more pronounced maximum at g™ than the experimental curves /(g), as already noticed 
in the literature 0, |23|- There are obviously large scale inhomogeneities in actual aerogels that are not reproduced 
in the simulations. Moreover, as emphasized in Refs.fT^ [2^. there are some significant details that are neglected in 
the DLCA model, such as the rotational diffusion of the aggregates, their polydispersity and irregular form, and all 
kinds of possible restructurating effects. 



B. Fluid structure during adsorption 

As shown in Refs.P, 0, the elementary condensation events (avalanches) that occur in 87% porosity aerogel as 
the chemical potential is slowly varied are always of microscopic size, whatever the temperature. This implies that 
the adsorption isotherms are smooth in the thermodynamic limit or when averaging over a large number of finite 
samples, as illustrated in Fig. 2. We have computed the correlation functions and the corresponding structure factors 
for a number of points along the T* — 0.5 and T* = 0.8 isotherms, as indicated in the figure. We first consider the 
lowest temperature. 

Figs. 3 and 4 show the evolution of the correlation functions hff{r) and hgf{r) with chemical potential. One can 
see that the curves change significantly as /i increases (note that the vertical scale in Fig. 3(a) is expanded so as 
to emphasize the presence of a shallow minimum in the curves; accordingly, the values of hff{r) near zero are not 
visible). For very low values of the chemical potential (e.g. fi = —6), hff{r) looks very much like hgg{r), apart from 
a shift towards larger values of r by a distance of about 2 lattice spacings. Indeed, as shown in our earlier workP,!^], 
in the early stage of the adsorption process, the adsorbed fluid forms a liquid fllm that coats the aerogel strands 
and whose thickness is approximately one lattice spacing at low temperature. In consequence, the distribution of the 
fluid particles follows the spatial arrangement of the aerogel, a feature already observed in a continuum version of 
the model|23. The existence of the liquid film also reflects in the rapid decrease of hgf{r) with r, which indicates 
that the fluid is only present in the vicinity of the gel particles (the fact that hgf (d) < hgf (a) may be ascribed to 
the connectivity of the gel: around a gel particle, there are always other gel particles - 2.5 in average in the first shell 
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FIG. 2: Average hysteresis loops in a 87% porosity aerogel at T* = 0.5 and 0.8 (from left to right). The points along the 
adsorption and desorption isotherms at T* = 0.5 indicate the values of the chemical potential for which the correlation functions 
are computed. The desorption isotherm has been computed either in presence of an external reservoir (solid line) or by using 
the procedure described in section IIIC (dashed line) (Color on line). 



- and the probability to find a fluid particle is thus suppressed). As /i increases, the magnitude of the fluid-fluid 




FIG. 3: Fluid- fluid correlation function hff{r) along the adsorption isotherm in a 87% porosity aerogel at T* — 0.5. (a): From 
top to bottom, the curves correspond to points 1 to 8 in Fig. 2; the dashed line is the gel correlation function hgg{r). (b): 
Magnification of (a) showing the evolution of the minimum as /i varies from —4.55 to —4 (points 5 to 8 in Fig. 2) (Color on 
line). 

correlations decreases at small r (the contact value decreases from 4.45 to 0.15) and the depth of the minimum in 
hff{r) decreases as it shifts to larger values of r (its location varies from 5 to 24 as /i increases from —6 to —4.47). The 
minimum disappears as the last voids in the gel fill with liquid (note the difference in the vertical scales of Figs. 3(a) 
and 3(b)), and finally, as one reaches saturation (at fXsat = —4), the shape of the gel-gel correlation function hggir) 
is approximately recovered, in accordance with Babinet principle |26ll29| . A similar behavior is observed in the cross 
correlation function hgf(r) in Fig. 4, but the minimum occurs at a smaller distance (its location is approximately 
r « 12 when it disappears). 

Like for ggg(r), we may associate to the location of the mimimum in hff(r) a length that characterizes the 
correlations within the adsorbed fiuid. The fact that becomes significantly larger than as the adsorption proceeds 
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shows that the fluid develops its own complex structure that does not reflect anymore the underlying gel structure. 
This is likely in relation with the fact that some of the condensation events studied in our previous works HH extend 
much beyond the largest voids in the aerogel^3- It is worth noting that these large avalanches (with a radius of 
gyration Rg « 12 2, 3]) occur approximately in the same range of chemical potential (—4.5 < fi < —4.4) where ^/ 
reaches its maximum (this also corresponds to the steepest portion of the isotherm). 




FIG. 5: Fluid-fluid structure factor Sff(q) along the adsorption isotherm in a 87% porosity aerogel at T* — 0.5. The numbers 
refer to points 1 to 8 in Fig. 2. Tlie dashed line is the gel structure factor Sgg{q) and the dotted line illustrates the influence 
of the continuation procedure for q <2n/L (see Jilj) (Color on line). 

The corresponding structure factors Sff{q) and Sgf{q) are shown in Figs. 5 and 6, respectively (sj. The main 
feature in Sff{q) is the presence of a broad peak that grows and moves towards smaller wavevector as the fluid 
condenses in the porous space. This peak is clearly associated to the minimum in hff[r) (its location is approximately 
proportional to ^) ^^^d it thus tells us the same story: the growing of a characteristic length scale in the fluid along 
the capillary rise. The peak disappears in the last stage of the adsorption process and is then replaced by a plateau 
(see curve 7 in Fig. 5). Finally, at ^ = Msat, one recovers a structure factor that can be deduced from Sgg{q) — 1 by 
a linear transformation, in accordance with Babinet principle (there are no oscillations in Sff{q), however, because 
the fluid- fluid hard-core diameter is zero). 

The evolution of the gel- fluid cross structure factor Sgf{q) is somewhat different. The peak is more prominent, 
as a consequence of the 'no- fluctuation' condition Sgf{q = 0) = (this feature may not be so pronounced in actual 
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FIG. 6: Same as Fig. 5 for the cross structure factor Sgf{q). Note that the vertical scale is not logarithmic because Sgf{q) has 
negative values (Color on line). 

systems because of large-scale fluctuations), and it is located at a larger wavevector than in Sff{q) (in line with the 
corresponding locations of the minima in hff{r) and hgf{r)). The most substantial difference with Sff{q) is that 
the amplitude of the peak starts to decrease much before the end of the adsorption process. The negative correlation 
observed at saturation is again due to complementarity with the gel-gel structure 13J. 

We have repeated these calculations at T* — 0.8 in order to investigate the influence of temperature. In the 87% gel, 
T* = 0.8 is just below Th, the temperature at which the hysteresis loop disappears (see Fig. 2). hff{r) and hgf{r) still 
exhibit a minimum that moves towards larger r upon adsorption. However, the characteristic length associated 
to the minimum of hff{r), does not exceed 14 lattice spacings, indicating that the size of the inhomogeneities in the 
fluid decreases with increasing T. A similar observation was made in Ref.0,|3| concerning the size of the avalanches 
which become more compact at higher temperature and often correspond to a condensation event occuring in a single 
cavity of the aerogel. The shape of the corresponding structure factors does not change significantly with respect to 
the T* = 0.5 case, but the amplitude is significanly reduced: the maximal amplitudes of the peaks in Sff{q) and 
Sgf{q) are divided approximately by 5 and 2, respectively. 

As shown in Refs. 0,0,13, temperature has a much more dramatic influence on the adsorption process in gels of higher 
porosity. In particular, at low enough temperature (T < Tc{(j)) with T*{(j)) « 0.5 in the 95% gel|3|), a macroscopic 
avalanche occurs at a certain value of the chemical potential, with the whole sample filling abruptly, which results in 
a discontinuous isotherm in the thermodynamic limit. In a finite system, the signature of a macroscopic avalanche 
is a large jump in the fluid density whose location in ^ fluctuates from sample to sample, which results in a steep 
but smooth isotherm after an average over the gel realizations (one then has to perform a finite-size scaling study to 
conclude on the proper behavior in the thermodynamic limit P,0|). Within a grand canonical calculation, there is 
unfortunately no way to study the evolution of the structural properties of the fluid during a macroscopic avalanche as 
this would require to consider intermediate fluid densities that are inaccessible ^2] (the situation would be different if 
the fluid density was controlled instead of the chemical potential, as is done frequently in experiments (^,33]). All we 
can do is to study the 95% gel at a higher temperature where the adsorption is still gradual, for instance at T* = 0.8. 
In this case, no qualitative change is found with respect to the case of the 87% gel at T* = 0.5. Indeed, as emphasized 
in Ref.Q, adsorption proceeds similarly in a high-porosity gel at high temperature and in a lower-porosity gel at low 
temperature. The correlation length is somewhat larger in the 95% gel (beyond 30 lattice spacings) so that finite- 
size effects come into play (in particular, it becomes problematic to extraplotate Sff{q) to g = so as to simulate 
the infinite-size limit). To go to lower temperatures, it would be thus necessary to use a much larger simulation box, 
which would increase considerably the computational work. Note that one expects hff(r) to decay algebraically at 
the critical temperature T(.{4>). Indeed, according to the analogy with the T = nonequilibrium random-field Ising 
model (RFIM), there should be only one important length scale in the system close to criticality, length scale which is 
proportional to the average linear extent of the largest avalanches [s^. At criticality, this correlation length diverges. 
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C. Fluid structure during desorption 

As discussed in detail in Refs. Jj,i^J, different mechanisms may be responsible for gas desorption in aerogels, depend- 
ing on porosity and temperature. In contrast with adsorption, the outer surface of the material where the adsorbed 
fluid is in contact with the external vapor may play an essential role. For instance, in the 87% aerogel at T* = 0.5, 
the theory predicts a phenomenon akin to percolation invasion : as /i is decreased from saturation, some gas "fingers" 
enter the sample and grow until they percolate at a certain value of fi, forming a fractal, isotropic cluster. The des- 
orption then proceeds gradually via the growth of the gaseous domain. Accordingly, in the thermodynamic limit, the 
isotherm shows a cusp at the percolation threshold followed by a steep but continuous decrease (the cusp is rounded 
in finite samples). 

The simulation of the desorption process thus requires the use of an explicit external reservoir adjacent to the 
gel sample, which of course introduces a severe anisotropy in the model and makes it difficult to calculate radially 
averaged correlation functions. To circumvent this problem, we have used another procedure where the desorption 
is not initiated by the interface with an external reservoir but triggered by the presence of gas bubbles inside the 
material. We have indeed observed in our previous studies (see Fig. 16 in Ref.Q) that the last desorption scanning 
curves (obtained by stopping the adsorption just before saturation and then decreasing the chemical potential) look 
very much like the desorption isotherms obtained in presence of a reservoir (when averaging the fluid density deep 
inside the aerogel, which gives a good estimate of the isotherm in the thermodynamic limit). Near the end of the 
adsorption process, the remaining gaseous domain is composed of isolated bubbles which obviously play the same role 
as an external reservoir when the chemical potential is decreased. The advantage of initiating the desorption with 
these last bubbles is that one can use the same geometry as during adsorption, with periodic boundary conditions in 
all directions (moreover, using small bubbles instead of a planar interface of size considerably suppresses finite-size 
effects). 

In practice, the calculation has been performed by keeping five bubbles in each sample. This implies that the 
chemical potential at which desorption is initiated is slightly different in each sample (if one chooses the same fj, 
in all samples, some of them may be already completely filled with liquid). This number of bubbles results from a 
compromise: on the one hand, keeping a single bubble may not be sufficient to trigger the desorption process (in some 
samples, the growth of the bubble is hindered by the neighboring gel particles and the desorption occurs at a much 
lower value of the chemical potential); on the other hand, keeping too many bubbles results in a too rapid growth 
of the gas domain. As can be seen in Fig. 2, the isotherm obtained with this procedure is indeed very close to the 
isotherm calculated in presence of an external reservoir. 

The fluid-fluid and solid-fluid correlation functions computed at several points along the desorption branch are 
shown in Figs. 7 and 8. 
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FIG. 7: Fluid-fluid correlation function hff{r) along the desorption isotherm in a 87% porosity aerogel at T* = 0.5. The 
numbers refer to the points in Fig. 2; the dashed line is the gel correlation function hgg{r). The vertical scale is expanded in 
the inset, showing the very slow decrease of hff{r) towards zero in the steepest portion of the isotherm (Color on line). 

One can see that hff{r) is dramatically changed with respect to adsorption: from saturation down to ^ « —4.7 
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FIG. 8: Same as Fig. 7 for the cross correlation function hgf(r) (Color on line). 

(curves 9 to 16), the function is monotonically decreasing, showing no minimum. Moreover, as shown in the inset 
of Fig. 7, a long-range tail is growing and hff{r) may differ significantly from zero at r = _L/2'36]. Although it is 
difficult to characterize this decrease by a unique and well-defined correlation length (the curve cannot be fitted by a 
simple function such as an exponential), it is clear that the range of the correlations is small at the beginning of the 
desorption process, then increases considerably, goes through a maximum in the steeppest portion of the isotherm 
(corresponding approximately to point 13 in Fig. 2), and eventually decreases. As the hysteresis loop closes and the 
adsorbed phase consists again of a liquid film coating the areogel strands, a shallow minimum reappears in the curve, 
which is reminiscent of the underlying gel structure. 

In contrast, the changes in the cross-correlation function hgf(r) between adsorption and desorption are very small 
(the function has only a slightly longer tail during desorption). It appears that the gel-fluid correlations depend 
essentially on the average fluid density: for a given value of p/, they are almost the same on the two branches of the 
hysteresis loop. 

The calculation of the fluid-fluid structure factor Sff{q) is complicated by the fact that hff{r) decreases very slowly 
to zero, and one cannot use anymore the continuation procedure that forces the sum-rule, Eq. 11, to be satisfied (the 
results for q < 2ti/ L change considerably with L, which shows that the resulting curve is not a good approximation of 
the infinite-size limit). Instead, we have only subtracted from hff{r) its value at r = L/2 (setting hff = hff{r = L/2) 
in Eq. 12) so as to avoid the large but spurious oscillations in Sff{q) that result from the discontinuity at L/2 (there 
is still a discontinuity in the slope of hff (r) at L/2 that induces small oscillations in some of the curves of Fig. 9). It 
is clear that flnite-size effects are important in this case and the results for q < 2t:/L must be considered with caution. 

The resulting fluid-fluid structure factors along the desorption isotherm are shown in Fig. 9. (We do not present 
the curves for Sgf{q) as they look very much like those in Fig. 6 with only a slightly broader peak.) As could be 
expected, the structure factors computed just before and after the knee in the isotherm (for /i > —4.67) are very 
different from those obtained during adsorption. Firstly, the peak that was associated to the minimum in hff (r) has 
now disappeared and the small-g intensity saturates to a value that is considerably larger than the maximum value 
obtained during adsorption (compare the vertical scales in Figs. 5 and 9). Secondly, as /i varies from —4.65 to —4.67 
(curves 11 to 14), there is a linear portion in Sff{q) whose maximal extension is about one decade on a log- log scale. 
On the other hand, when /i is decreased further, the peak in Sff{q) is recovered and the curves become more similar 
to the ones computed on the adsorption branch. 

The linear regime in Sff{q) strongly suggests the presence of fractal correlations |35j . However, according to our 
previous studies 2], it is only the gaseous domain that should be a fractal object at the percolation threshold, as 
illustrated by the isotropic and strongly ramifled structure shown in Fig. 10. The correlation function hff{r), on the 
other hand, does not discriminate between a site representing a gaseous region [pi « at low temperature) and a 
site representing a gel particle {pi = 0). In order to really show the existence of fractal correlations within the gas 
domain, one must consider either the 'gas-gas' correlation function (deflning the quantity pf"** = rji — pi which is equal 
to 1 only in the gas phase) or the complementary function that measures the correlations within the dense (solid or 
liquid) phase (one then deflnes pj'^^'^'^ = l — pf^^). The corresponding structure factor Sdd{q) is the quantity that is 
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FIG. 9; Fluid- fluid structure factor Sff{q) along the desorption isotherm in a 87% porosity aerogel at T* — 0.5. The numbers 
and arrows refer to points 9 to 18 in Fig. 2. The dashed line is the gel structure factor Sgg{q) (Color on line). 




FIG. 10: Snapshot of the vapor domain in a 87% gel sample during desorption at T* — 0.5 and /i = —4.63 (Color on line). 



measured experimentally when using the 'contrast matching' technique|^. It is related to Sff{q) and Sgf{q) by 

(Pg + Pf){Sdd{q) - 1) = Pf{Sff{q) - 1) + ^PgPfSgf{q) + Pg{Sgg{q) - 1) . (15) 

Sdd{q) is shown in Fig. 11 for ^ — —4.65 37], in the region of the knee in the desorption isotherm (point 11 in Fig. 
2). It clearly contains a linear portion over almost one decade and can be very well represented in this range of 
wavevectors by the fitjs^ 



Sdd{q) 



sin [{df - l)tan-i(gZ)] 



2l('i/-l)/2 



(16) 



with df = 2.45 and I = 17, where Z is a crossover length that limits the fractal regime at large distances. Note that the 
linear portion itself has a slope « — 2.1 (the above formula reproduces the right slope only when I is very large p^). 
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An accurate determination of would therefore require a much larger system and, at this stage, it is not possible 
to decide if the fractal dimension is consistent with that of random percolation. In any case, these results strongly 
suggest that the gas domain exhibits fractal correlations during desorption, correlations which have no relation with 
the underlying gel microstructure (we recall that there is almost no fractal regime in the 87% aerogel, as can be seen 
in Fig. 1). 
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FIG. 11; Structure factor of the dense phase (see text), Sdd{<l), during desorption at T* = 0.5 and fj, = —4.65. The dashed-dotted 
curve is the best fit according to Eq. 16 and the straight dashed line has a slope —2.1. 

Raising the temperature to T* = 0.8 has a dramatic effect, as shown in Fig. 12. The maximum value of Sff{q) has 
dropped by two orders of magnitude and there is no significant region with a fractal-like power-law behavior. Indeed, 
hff{r) has no more a long-range tail and the correlations are very similar during adsorption and desorption, as could 
be expected from the very thin shape of hysteresis loop. This is the signature that the desorption mechanism has 
changed, in agreement with the analysis of Refs.'l','5|. It is now due to a cavitation phenomenon in which gas bubbles 
first appear in the largest cavities of the gel and then grow and coalesce until the whole void space is invaded |39j. 
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FIG. 12: Fluid-fluid structure factor Sff(q) along the desorption isotherm in a 87% porosity aerogel at T* = 0.8. (Color on 
line) . 

We have not studied the correlations during desorption in the 95% porosity gel. At very high temperature (T* > 0.9), 
desorption is expected to be due again to cavitation, 2], and the results shoud be similar to those in the 87% that 
have just been described. On the other hand, at low temperature (e.g. T* = 0.5), the theory predicts a depinning 
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transition in which a self-affine interface sweeps through the whole sample, resulting in a discontinuous desorption 
isotherm^]. Therefore, like in the case of the macroscopic avalanche during adsorption, the correlations along the 
isotherm cannot be studied within the framework of a grand-canonical calculation. At intermediate temperatures, one 
could probably observe again extended fractal correlations associated to a percolating cluster of gas, but this study 
requires the use of larger systems so as to probe smaller values of q and discriminate the effects due to the own fractal 
structure of the gel. 

IV. SCATTERED INTENSITY AND COMPARISON WITH EXPERIMENTS 

As mentionned in the introduction, there have been two recent scattering studies of gas condensation in aerogels, 
both with "'He^, 1^. In Ref.Q, light scattering is used to study adsorption and desorption in a 95% porosity gel 
at several temperatures between 4A7K {T* « 0.86) and 5.08K [T* k, 0.98). These experimental results cannot 
be directly compared to our theoretical predictions, and our system size is too small to investigate the large-scale 
inhomogeneities that are seen in the experiments (some of them are visible to the eye). However, there are two key 
observations that appear to be in agreement with our predictions: i) at the lowest temperature studied, the optical 
signal due to helium adsorption is larger than if the fluid density was simply correlated to the density of silica, 
indicating that the correlations within the fluid extend beyond the aerogel correlation length, and ii) the aerogel is 
much brighter during desorption, indicating that the characteristic size of the density fluctuations is much larger than 
during adsorption. 

This latter conclusion was also reached from the small-angle x-ray scattering measurements (SAXS) performed 
in a 98% porosity aerogel at 3.5if [T* k, 0.67) j^]. SAXS is particularly well suited for observing the structural 
features associated to fluid adsorption, and in order to compare more easily to experiments we shall push further our 
calculations and compute the resulting scattered intensity. Of course, the predictions must be taken with a grain 
of salt, considering the limitations in the model and the theory. Since the scattered intensity is proportional to the 
Fourier transform of the electron density fluctuations, one has 

I{q) (X pgF{qfSM + 2^p;pJaFiq)Sgf{q) + Pfa^Sff{q) (17) 

where F{q) is the form factor of silica particles (see section IIIA) and a is the ratio of the electron density in the 
adsorbed fluid to that in the solid. As an additional approximation, we shall take F{q) « 1, restricting the study 
to the range 2'k/L < q < 2 where this is presumably a reasonable approximation (in real units this corresponds to 
0.02 < g < 0.7 nm~^, taking a — 3nm). Assuming that the adsorbed hquid has the same density as the bulk hquid 
at Psat; and using the tabulated densities of helium and silica, one finds that a varies from 6.55 10~^ at T* = 0.5 to 
5.75 10-2 at T* = 0.8. 

The theoretical scattered intensities during adsorption and desorption in the 87% gel at T* = 0.5 are shown in 
Figs. 13 and 14. Like in the experiments^, we plot the ratio R{q) = I{q)/Ie{q) where Ie{q) is the contribution of the 
empty aerogel (the first term in the right-hand side of Eq. 17) in order to accentuate the effects of the adsorbed gas, 
especially in the initial stage of adsorption. Moreover, we hope that this also partially corrects the small-(? defects 
due to the absence of large-scale fluctuations in the DLCA gel structure factor. 

The main features of the curves displayed in Fig. 13 are the following: i) At the very beginning of adsorption, R{q) 
slightly increases but remains almost independent of q. This is the signature of the ^He film coating the aerogel. In 
this regime, the main contribution to the scattered intensity comes from the gel-fluid correlations (the second term in 
the right-hand side of Eq. 17). ii) As p, increases, the scattering grows in intensity at small q, reflecting the presence 
of the broad peak in Sff{q) that moves towards smaller wavevector with filling (see Fig. 5). iii) As the aerogel fills 
further, R{q) decreases until it becomes again almost flat at complete filling. The total intensity is then reduced with 
respect to that of the empty aerogel. Direct comparison with the experimental results of Ref.jg] is again problematic: 
the adsorption isotherm in the 98% gel is indeed very steep at 3.5K, suggesting that one may be in the regime of a 
macroscopic avalanche. However, the behavior of the experimental R{q) is remarkably similar to what have just been 
described (preliminary measurements in a 86% aerogel l405 also show the same trends). The results of Ref.Jij were 
interpreted according to a model of two-phase coexistence, with a 'film' phase in equilibrium with a filled 'pore' phase. 
This is at odds with the theoretical scenario discussed in Refs.P,|3,l3| which emphasizes the nonequilibrium character 
of the transition. The present results seem to show that this approach can also elucidate (at least qualitatively) the 
behavior of the scattered intensity. 

During desorption, the most characteristic feature of the curves shown in Fig. 14 is the very significant increase of 
the ratio R{q) at small q with respect to adsorption (note the logarithmic scale on the vertical axis). This is related 
to the corresponding increase in Sff{q) shown in Fig. 9 and is clearly due to the presence of long-range correlations 
within the fluid. As the desorption proceeds, R{q) goes through a maximum and then decreases until it becomes flat 
again. Remarkably, no power-law fractal regime is visible in R{q) in the range 0.06 ^ 9 ^ 1 as was the case with 
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FIG. 13: Theoretical ratio R{q) of the scattered intensity 7(g) to the scattered intensity leil) of the empty aerogel during 
helium adsorption in a 87% aerogel at T* = 0.5. The dashed line coresponding to f? = 1 is shown for reference (Color on line). 
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FIG. 14: Same as Fig. 13 during desorption (Color on line). 

Sdd{q) in Fig. 11. It is the small value of a (due to the small electron density of He), and not the division by Ie{q), 
which is responsible for this unfortunate disappearance (/(g) becomes proportional to Sdciil) when a — 1, which is 
only the case in a contrast matching experiment). In the measurements of Ref.|^, this increase of R{q) at small q is 
not mentionned, but the analysis of the data shows that the charateristic size of the inhomogeneities is much larger 
than during adsorption, as already mentionned, and that it decreases rapidly in the last stage of desorption. 

Not surprisingly, the theoretical scattered intensity in the 87% aerogel is considerably smaller at high temperature, 
as illustrated in Fig. 15 for T* — 0.8. The intensity ratio R{q) has been divided by about 40. We therefore conclude 
that the magnitude of the scattered intensity can indicate that the nature of the desorption process has changed. We 
leave to our colleague experimentalists the challenge of checking the presence of a fractal-regime during desorption, 
as was done in Vvcor |10nlli.il2j and xerogel|14j. 
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FIG. 15: Same as Fig. 13 during desorption at T* = 0.8 (Color on line). 
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